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ABSTRACT 

Wc present the evolutionary models of metal-free stars in the mass range from 0.8 to 1.2M Q 
with up-to-date input physics. The evolution is followed to the onset of hydrogen mixing into a 
convection, driven by the helium flash at red giant or asymptotic giant branch phase. The models 
of mass M > 0.9Mq undergo the central hydrogen flash, triggered by the carbon production due 
to the 3a reaction. We find that the border of the off-center and central ignition of helium core 
flash falls between 1.1 and 1.2 Af Q ; the models of mass M < 1.1M Q experience the hydrogen 
mixing at the tip of red giant branch while the models of M = 1.2M© during the helium shell 
flashes on the asymptotic giant branch. The equation of state for the Coulomb liquid region, 
where electron conduction and radiation compete, is shown to be important since it affects the 
thermal state in the helium core and influences the red giant branch evolution. It is also found 
that the non-resonant term of 3a reaction plays an important role, although it has negligible effect 
in the evolution of stars of younger populations. We compare our models with the computations 
by several other sets of authors, to confirm the good agreement except for one study which 
finds the helium ignition much closer to the center with consequences important for subsequent 
evolution. 

Subject headings: stars: general — stars, Pop III, evolution, red giant 



1. Introduction 

Although the initial mass function of metal- 
free stars formed out of primordial matter has 
not yet been determined, there is evidence that 
metal-free stars of low and intermediate mass were 
formed. For example, the fact that the frequency 
of stars with strong carbon enhancements is much 
larger than in the case of population I an d II 
stars ( Rossi et "ail 1 1998 ; Lucatello et al.ll2005l see 
also the review by Beers & Christlieb 2005) is 
consistent with theoretical predictions. One of 
the most prominent characteristics of models of 



low-mass and intermediate-mass extremely metal- 
poor stars (EMPSs) is that they become car- 
bon stars at an earlier stage of evolution than 



the s tars of younger populations (jFuiimoto et al 
2000l hereinafter FII00). Recently, it is argued 
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that EMP stars survived to date were formed as 
a secondary compon ent in the binary sy stems of 
massive companions (jKomiva et al.ll20071 ). As an- 
other example, the peculiar abundance character- 
i stics of the most met al poor stars yet discovered 
( Christlieb et al .120021) can be understood in terms 
of the evolutionary properties of Z = or ex- 
tremely iron-poor models in an interm ediate-mass 
interacting binary (jSuda et alj 120041 ). Thus, the 
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evolutionary characteristics of Z = models have 
direct relevance to discussions of star formation in 
the early universe. 

EMP stars may alter their initial surface abun- 
dances by bringing to the surface products of inter- 
nal nuclear transformations which occur in unique 
ways. For initial CNO abundances Zqno S 10~ 7 , 
the outer edge of the convective zone, generated 
by the first off-center helium flash in the hydrogen- 
exhausted core extends into the hydrogen-rich lay- 
ers, eventually leading to the enrichment of carbon 
and nitrogen in the surface (Fujimoto et al. 1990, 
hereinafter FIH90; Hollowell et al. 1990). This 
"He-flash driven deep mixing (He-FDDM)" mech- 
anism is distinguished from the third dredge-up 
in asymptotic giant branch (AGB) stars of Popu- 
lations I and II wh ich enriches surface material 
in 12 C (jlbenl I1975T ). Even if we take into ac- 
count the possible surface pollution by the accre- 
tion of metal-rich gas after birth, the interior evo- 
lution of metal-poor stellar models is not altered 
as long as the original CNO abund ance satisfies 
Zcno < 10" 8 (|Fuiimoto et al.lll995l ). 

Many calculations of Population III (here- 
inafter Pop. Ill) sta r evolution have been pub- 
lished (see FIH90; ISuda et all |2004 and refer- 
ences therein). Based on a top heavy initial mass 
functi on for the primordial cloud Ezer fc Cameron! 
()197ll ) first computed sets of massive metal-free 
models. Computations of the evolution of low 
ma ss metal - deficie nt stars were first carried out 
by IWagner (1974) in order to provide a set of 
stellar lifetimes for use in calculations of galactic 
chemical evolution. Computations of low-mass, 
zer o-metallicity star evol u tion were performed 
by ICastellani fc Paolicchi ( 1975 ) up to the ex 



hausti on of hydrogen at the center. iD'Antona 
(|1982l ) pursued the evolution of metal-free mod- 
els of mass 1 Mq to the red giant branch (RGB) 
phase, but because of a large time step in her 
co mputation, did not find any a bnormal bchav 



Guenther fc Demarque ( 1983[ ) suggested that 



a convective instability occurs during the cen- 
tral hydrogen-burning phase and FIH90 found 
this to be the case. The "peculiar" evolution 
of Z = model stars and its interpretation 
is revealed by the c alculations of FIH90 and 



Hollowell et al 



( 19901) . who followed the evolu- 
tion of a Z = star of mass 1M Q and gave 
explanations for the core He-H flash, a shell 



He-H flash, and the He-FDDM phenomenon. 
By computing in detail the progress of the 
hydrog en-mixing event and t he subsequent evo- 
lution, Hollowell et al. ( 1990[ ) first demonstrated 
that metal-free stars b ecome CN rich carbon stars 
at t he red giant stage. Cassisi fc Castellani ( 19931 ) 
and lCassisi et al.l (|1996l ) found the core He-H flash 
in models of Z = 10~ 10 , but, since their compu- 
tations did not continue beyond the initiation of 
the helium core flash, they found no evidence for 
the ad ditional events just described. IWeiss et al 



( 2000f ) (hereinafter WOO) also made computations 
of low mass metal-free stars, but terminated their 
computations just after the ignition of the helium 
core flash while the flash-driven convective zone 
was still growing in mass. 

hcrc- 



Computations by ISchlattl et alj (|2001 



inafter SOI) support the reality of the He-FDDM 
phenomenon for low-ma ss zero metallicity stars, 
where as computations bv lSiess. Livio fc Lattanziol 
(|2002t ) (hereinafter SLL02) do not. While they 
found the mixing of hydrogen into the helium- flash 
driven convective zone, the failure of SLL02 to re- 
produce the He-FDDM phenomenon may be as- 
cribed to their assumption of instantaneous mix- 
ing of elements in the helium flash convective zone; 
in consequence, nuclear energy released from the 
mixed hydrogen is distributed throughout the en- 
tire convective zone and the entropy of the zone is 
not built up sufficiently for the edge of the convec- 
tive zone to reach deeply into the hydrogen profile. 
Because mixing and burning of hydrogen occur on 
a very short timescale in the middle of the helium 
convective zone, following the He-FDDM event 
properly requires a trea t ment of time-dependent 
mixin g (jHollowell et al.l Il990h . iHollowell etTaf 



(1990) lirst derived the surface composition caused 
by the He-FDDM p henomenon for a Z = model 
Schlattl et all (|2002h also de- 



o- 



of mass M = 1M ( 
rived the surface composition for a model of mass 
M = 0.82M Q and compared the results with the 
observations of EMPSs. 



Aft er the discovery of HE0107-5240 (jChristlieb et al 
20021 ). which held the distinction of being the 



most metal-poor gia nt known before th e discov 



cry of HE1327-2326 (|Frebel et al.l l2005h . models 



of EMPSs were computed in an effort to deter- 
mine the evolutionary h istory of the observed 



star Both iPicardi et all (|2004 hereinafter P04) 



and lWeiss et alj (|2004l ) found the He-FDDM phe- 
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nomcnon for Z = models of mass M = 0.8Mg 
and M = 0.82M Q , respectively. 

All but one of the works just cited find the 
core Hc-H flash; thus far, only Fujimoto and his 
coworkers have found the shell Hc-H flash driven 
by a violent CN-cycle burning at the base of the 
hydrogen-burning shell. Why various groups find 
different evolutionary characteristics for similar 
masses and compositions has not yet been fully 
understood. One reason for this is that the com- 
putation of Z = models requires not only a care- 
ful treatment of the input physics but also a care- 
ful numerical solution. In particular, the treat- 
ment of convective regions requires consideration 
of many complicated factors including uncertain 
parameters. SOI made comparisons among mod- 
els by changing parameters related to abundances, 
diffusion, stellar mass, and the mixing length al- 
gorithm for convection, and suggestions have been 
made for the causes of differences such as radiative 
and/or conductive opacities and neutrino energy- 
loss rates (WOO; SOI; SLL02). However, we are 
not satisfied that the basic causes for differences 
have been pinpointed satisfactorily. 

In this paper, we examine the differences in the 
evolutionary characteristics of low-mass Z = 
models by taking into account the choice not only 
of the radiative and conductive opacity and of neu- 
trino energy-loss rates, but also of nuclear reaction 
rates. Thanks to the many efforts until today, use- 
ful sets of input physics are available and they 
have been more and more accurate for numeri- 
cal computations. But we think it is worth try- 
ing to check the dependence on the input physics 
and to elucidate the structural characteristics in 
detail because these efforts are still going on and 
because the evolution of low-mass Popualtion III 
stars has not been established yet. In particular, 
as for the three alpha reaction rates, the determi- 
nation of resonance states in 12 C is still controver- 



sial (jFvnbo et al.l 120051 ); the simulations of three 
body interactions are challenging, and yet, become 
feasible task in near future tha nks to the improved 
comp uter resources (see e.g., iKurokawa fe Kat^" * 



20051 . and the references there). 

In the next section, we elaborate the input 
physics adopted in the stellar evolution program. 
In §3, we present the result of computations of evo- 
lution of low mass Z = model stars with updated 
input physics, and discuss the model characteris- 



tics and their dependences on the input physics 
including the resonant and non-resonant reaction 
rates of 3 a reactions. In section 4, comparisons 
with other works with the different input physics 
taken into account, and, in §5, we summarize con- 
clusions. 

2. The Computational Program 

The original program to com p ute st ellar evo- 
lution was constructed by llben ( 1965|) and has 
been modified p eriodically (e.g., see llben 1975: 
Iben et al.lll992f ). The size of each mass shell and 
the time step are regulated, respectively, by the 
gradients with respect to space and time of struc- 
ture and composition variables. Typically, 200-300 
mesh points are required for main sequence models 
and 700-800 for core helium-burning models with 
hydrogen-burning shells. The number of equilib- 
rium models required to follow evolution from the 
zero-age main sequence to the end of RGB or AGB 
phase varies from 5000 and 30000, with the ex- 
act number depending on initial mass. The stellar 
structure equations are typically satisfied to better 
than one part in 10 5 . 

In the program, the radiative opa city is ob- 

tained by interpolation in OPAL tables (llglesias fc Rogers! 
19961) and in tables by Alexander fe Ferguson! 
1994 hereinafter, AF 94 tables) and th e con- 



ductive opacity is from Itoh et all (1983f) . The 



equation o f state invo l ves fits bv | Iben et al.l (1199 



to work bvlAbel dl959h iBowers fc Salpeterl (U96 




Slattery et all (Il980h.lsiatterv et all (Il982h llvetomi fc Ichimaru 
(ll982l).lHanserJ (ll973MH"ansen fc Maz!gh3 (Il978l). 
Cohen fc Kefferl (jl955l ). and ICarrl (Il96lh Neu- 
trino energy-loss rates are from lltoh et al.l ( 19961 
in the following 196) for plasma, photo, pair 
and bremsstrahlung processes. In order to com- 
pare with other works, we make use of two 
sets of nuclear reaction rates: those given by 
Caughlan fc Fowler! (|l988l . in the following CF88) 
and those given in th e latest NACRE compilation 



( Angulo et al. 19991 ). Nuclear screening factors 



for weak and strong screening are also taken into 
cons ideration using stan dard prescriptions (see, 
e.g.. lBohm-Vitenselll958l ). but only weak screen- 
ing is dominant in the actual computational range 
in this work. Nine nuclear species are considered: 
iH, 3 Hc, 4 He, 12 C, 14 N, 16 0, 18 0, 22 Nc, and 25 Mg. 
Abundances of these isotopes are determined by 
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16 nuclear reactions which include proton, elec- 
tron, and alpha captures. 

In regions of high temperature and high density 
which arc not covered by OPAL tables, we us e the 
analyt ical radiative opacit i es from the fits bv llbenl 
(|l975l ) to ICox fc Stewartl (|l970allbh ones. At ta- 
ble boundaries, no interpolation is made between 
table values and analytical values. Hence, jumps 
in the radiative opacity occur at table boundaries, 
but jumps are normally smaller than a factor of 2 
and do not seriously affect model convergence, pri- 
marily because, in regions not covered by the ta- 
bles, the overall opacity is dominated by electron 
conductivity At low temperature and low den- 
sity, interpolation between OPAL and AF94 tables 
is accomplished by setting k = kopal (1 — 0) + 
kaf 0, where = sin 2 (tt/2)(T - Ti)/(T 2 - T{) 
for the temperature T between T\ = 8000 K and 
T 2 = 10000 K, and k pal and k A f are OPAL 
and AF94 opacities, respectively. In surface layers 
of all models discussed here, AF94 tables provide 
opacities for all densities and temperatures en- 
countered. In regions of overlap, OPAL and AF94 
opacities agree rather well, so switching from one 
table to the other introduces little uncertainty. 

Although the quantum correction to the con- 
ductiv ity at low temper atures has been calcu- 



lated (jMitake et al.lll984h . it is now believed that 



the approximation employed is not appropriate at 
the temperatures considered. Therefore, we use 
the analytic fits devised by 183 for the conduc- 
tive opacity at low temperature; the liquid metal 
phase for various elemental compositions as well 
as strong degeneracy are taken into considera- 
tion. Strong degeneracy prevails if T <C Tp, where 
T F {K) = 5.930 xlO 9 ((l + l.OlS^e/^) 2 / 3 ) 1 / 2 -!), 
Pq is the density in units of 10 6 g cm -3 , and p e is 
the electron molecular weight. For the actual com- 
putations, since their approximations are good for 
T < Tp, we adopt the analytical approximations 



of lltoh et all 1)19831 ) for T < 0.5 T F . For T > 2 T F , 
we adopt th e conductive opaci ty used in the fits by 
Ibenl ( 1975 ) to calculations bv lHubbard fc Lampe 



1969 ) for non-rclativistic electrons and to lCanuto 
19701 ) for relativistic electrons. For intermediate 
temperatures, interpolation between the two ap- 
proximations is achieved by "sine square" weight- 
ing, analogous to that used in interpolating be- 
tween OPAL and AF tabular opacities. Finally, 
the conductive opacity is subjected to a limit de- 



scribed in the Appendix. Interpolated conductive 
opacities are given in Figure [IJ which covers the 
overall range of application for stellar evolution. 

For neutrino energy-loss rates, we adopt the 
fitting formulae presented by 196, who include 
pair, photo, plasma, and bremsstrahlung with 
weak degeneracy, the liquid-metal phase with low- 
tempcraturc quantum corrections and the crys- 
talline lattice phase, and the recombination neu- 
trino processes. Among these latter processes, we 
do not adopt the low-temperature co rrection to 
the liq uid-metal phase as calculated bv lltoh et al 



(|1984I ) since the corrections are generally small, as 
stated in 196. For T < 10 7 K, we assume that neu- 
trino energy-losses of all kinds can be neglected. 
The treatment of a multicomponcnt gas, which is 
important for bremsstrahlung, is considered in ap- 
plying the neutrino energy-loss rates of 196, as de- 
scribed in the Appendix. 

Mass loss is neglected in our work because we 
are interested primarily in the stellar interior and 
the evolution of the helium core is not affected by 
modest surface mass loss. To determine the tem- 



()l958l) ( 



see, e.g., 



length recipe by 


Bohm-Vitense 


Cox & Giuli 1968 


) and the mix- 



ing length is taken to be 1.5 times the local pres- 
sure scale height. Neither overshooting nor semi- 
convection are considered. 

In our calculations we focus on the evolution- 
ary trajectory up to the He-FDDM event and on 
the thermal structure of a star at stages of interest 
along its path to this event; our computations are 
terminated at the onset of hydrogen mixing into 
the convective zone driven by a helium flash. This 
event can occur during an off-center helium flash, 
and/or at the beginning of the thermally pulsing 
AGB (TPAGB) phase, depending on the initial 
mass and input physics. Subsequent evolution af- 
ter the mixing event is beyond the scope of this- 
paper and will be discussed in detail in a separate 
paper (Suda, Fujimoto, & Iben, in preparation). 

3. Evolution of Low-Mass Pop. Ill Stars 

The evolution of Z = stars has been com- 
puted from the zero-age main sequence (ZAMS) 
through the RGB and/or through the TPAGB 
phase for model masses of 0.8 — 1.2M©. For all 
models, the initial chemical composition is X = 
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0.767, Y = 0.233, and Z = 0, and the initial 
abundance by mass of 3 He is 2 x 10~ 5 . These 
abundances are based on models of big bang nu- 
cleosynthesis and are the same as those chosen by 
FII00. 

Table [T] lists the models computed with all the 
input physics updated for this paper. The first 
two columns give, respectively, the model identi- 
fier and the initial model mass; the labels "nac" 
and "cf" mean that the model has been com- 
puted with the choice of nuclear reaction rates 
of NACRE and CF88, respectively. The third to 
the ninth columns give the effective temperature 
and surface luminosity at the turn-off point and at 
the tip of the RGB, the times to reach these two 
stages, respectively. 

Figures [2] and [3] show, respectively, evolution- 
ary tracks in the H-R diagram and in the central- 
density and central-temperature plane for "cf" 
and "nac" models of mass in the range 0.8 — 
1.2Mq. All models are terminated at the on- 
set of hydrogen mixing episodes. We found that 
M < I.IM0 models ignite the helium burning in 
the off-center shell and the engulfment of hydro- 
gen by the flash-driven convection occurs just after 
the peak of the major core helium-flash at the tip 
of RGB. On the other hand, 1.2M Q models un- 
dergo the helium core flash at the center without 
the hydrogen mixing, and encounter the hydrogen 
mixing in the helium-flash convection during the 
helium shell flash at the beginning of thermally 
pulsating AGB (TP-AGB) phase. The two mod- 
els display almost the same trajectory irrespective 
of the choice of nuclear reactions, with the small 
differences stemming from the difference in the 3a 
reaction rates, less than a factor of 2 in the tem- 
perature range logT = 7.8 — 7.9. 

The evolutionary characteristics of the core he- 
lium flash and the hydrogen mixing event are sum- 
marized in Table [2] Each column gives the model 
identifier defined in Table [TJ the helium core mass 
Mf 18,1 and the helium burning rate Ljg ax when 
the helium burning rate reaches maximum, the 
mass coordinate Mg^g and the maximum tem- 
perature TgQg at the base of the convective shell, 
driven by the core helium flash, the helium burn- 
ing rate at the onset of hydrogen mixing at 
the RGB, and the time intervals, At' and Ai m j x , 
to it from the appearance of helium flash-driven 
convection and from the stage of maximum helium 



burning, respectively. The mass Mi of helium core 
is defined as the mass coordinate where the abun- 
dance of hydrogen is half of the surface abundance 
of hydrogen. 

In this section, we discuss in detail the evolu- 
tionary behavior of Pop. HI models with respect to 
the differences in initial mass and input physics. 

3.1. Hydrogen Burning Phase 

Variations in several quantities characterizing 
model stars are summarized in Figured] MS R func- 
tion of the hydrogen abundance at the center. In 
a Pop. Ill model, due to the absence of CNO cat- 
alysts, p-p chain reactions are initially the only 
mode of energy generation by hydrogen burning. 
Because of the weak temperature dependence of 
the energy-generation rate, the central tempera- 
ture keeps rising as the hydrogen abundance de- 
creases (top panel of Fig. [3}. As the central 
temperature increases, the 3a reaction gradually 
becomes active and CNO-cycle reactions begin to 
occur. Eventually, CNO-cycle reactions dominate 
the p-p chains with regard to total energy pro- 
duction. When this first occurs, the central abun- 
dance by mass Aqj^q of catalysts is I0~ n — I0~ 9 
at logT c ~ 7.8 — 7.9 (top and middle panels of 
Fig. HJ. Because of lower central temperatures, 
the less massive the star, the later is the evolu- 
tionary stage (and the smaller is the central hydro- 
gen abundance) at which the transition from burn- 
ing dominated by the p-p chains to burning domi- 
nated by the CNO-cycle reactions takes place. For 
M < 0.8M Q , hydrogen is depleted in the center 
before a transition can take place. After the CNO 
cycle takes over as the main source of energy gen- 
eration, because of the strong temperature depen- 
dence of CNO-cycle reactions, the central temper- 
ature remains nearly constant. At the same time, 
the core expands because of the central concen- 
tration of energy generation and the rate of pro- 
duction of catalysts slows down; the abundances 
of catalysts saturate at Aqj^q ~ I0~ 10 — I0~ 8 . 

The transition to the CNO-cycle dominated 
phase is accompanied by the formation of a con- 
vective zone which develops outward from the cen- 
ter as described in the lower panel of Fig. 2] The 
growth of the convective core is also evident in 
the middle panel of Fig. 21 which shows that the 
hydrogen abundance at the center, X c , stops de- 
creasing monotonically and increases for a time as 
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evolution progresses. Central convection is caused 
by a thermonuclear runaway and persists after the 
transition; due to the larger temperature depen- 
dence of the CNO-cycle energy-generation rate, 
energy generation is highly concentrated toward 
the center. Both the inward mixing of hydrogen 
from outer hydrogen-rich layers and the outward 
mixing of CNO elements generated near the cen- 
ter during the thermonuclear runaway amplify the 
average hydrogen-burning rate over the region en- 
compassed by convection relative to the average 
rate in the absence of convection. 

For stars of mass 0.9 < M/M Q < 1.2, the 
transition from the p-p chain dominated phase to 
the CNO-cycle dominated phase is delayed until 
the core has already begun to contract rapidly 
and electrons have begun to become degenerate 
at the center; the maximum energy generation 
has already shifted away from the center and 
core contraction has initiated the expansion of 
the envelope. A thermonuclear runaway called 
the core helium-hydro gen (He-H) flash takes place 
(|Fujimoto et al J Il990h ■ In Fig. [2 a first increase 
in the central temperature with decreasing cen- 
tral density indicates that the electron degeneracy 
is lifted; then, the temperature turns to decrease 
with the density so as to settle in the thermal equi- 
librium state where the nuclear energy generation 
balances the energy loss from the core. In this 
runaway, helium burning plays a key role through 
the production of CNO-cycle catalysts. Because 
of the strong temperature dependence of CNO- 
cycle reaction rates, the flash starts even when 
the contribution of the CNO cycle to the total 
energy-generation rate is smaller than that of p-p 
chain reactions in central regions. Smaller mass 
stars experience stronger flashes since the central 
entropy is smaller and the electron degeneracy is 
stronger at the onset of the thermonuclear run- 
away. Thus, the temperature reaches larger, lead- 
ing to a greater production of CNO catalysts and a 
greater extension of convection. For stars of mass 
M < 0.8M Q , the central temperature does not 
become large enough for the production of suffi- 
cient carbon to activate the CNO cycle, irrespec- 
tive of the choice of reaction rates. We see a small 
hump along the trajectories of M = O.8M mod- 
els in Fig. [3l which marks where the exhaustion of 
hydrogen at the center quenches a thermonuclear 
runaway. 



These evolutionary characteristics during core 
hydrogen burning are common to all of the models, 
regardless of the adopted nuclear reaction rates 
since the evolutionary tracks prior to the transi- 
tion are dominated by the p-p chain reactions with 
relatively small temperature dependences. All the 
models of masses in the range 0.9 < M/Mq < 1.2 
experience the core He-H flash. We may notice 
that the differences in the ignition temperatures 
are rather small because of the strong tempera- 
ture dependence of the 3a reaction rate, the flash 
grows weaker for models of larger mass. 

Core convection lasts until hydrogen is almost 
exhausted at the center (X c < 10~ 4 ). After the 
model settles in the thermal equilibrium, the cen- 
tral density increases again as the hydrogen abun- 
dance decreases. The central temperature at first 
rises, and then, decreases as the nuclear energy- 
generation rate decreases at the center (when 
X c < 2 — 3 x 10~ 4 ). Note that since the nuclear en- 
ergy generation is dominated by contribution from 
off-center burning, the model at this phase tends 
to have an isothermal core with the central tem- 
perature directly reflecting the temperature in the 
H-burning shell. In Fig. [3J we can see the rise 
in central temperature with the growth of core, a 
sudden jump in temperature occurring when the 
hydrogen burning shell passes the shells occupied 
by the flash convection, and hence, have larger 
CNO abundances. 

The resultant loop, as seen in Fig. [31 is char- 
acteristics of flash, implying that He-H core flash 
eventually exerts work through the expansion and 
contraction of core, during which the envelope first 
contacts and then expands again, as seen in Fig. [2] 
In the subsequent evolution, our models do not 
encounter the phenomena, so called "shell He-H 
flash" at the base of hydrogen burning phase, as 
FIH90 discovered. FIH90 discuss the behavior of 
hydrogen burning shell in the Z = environment 
and found the convective instability by the pro- 
duction of carbon and of nuclear energy by CN 
cycles. This is the counterpart of core He-H flash, 
i.e., electrons are degenerate at the base of hydro- 
gen burning shell, which causes He-H flash. How- 
ever, no other groups than Fujimoto and his col- 
laborators find such a event. The different conse- 
quence is due to the consideration of non-resonant 
effect in 3a reaction rate as we will discuss in 347 
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3.2. Helium Burning Phase 

Characteristic properties of the helium core 
flash are given in Tabled The thermal structure 
of the core when the helium core flash is ignited 
are determined by (1) neutrino energy losses which 
produce cooling, (2) the radiative and conductive 
opacity which controls heat flow, and (3) the tem- 
perature and energy generation in the hydrogen- 
burning shell which affects the heating due to 
gravitational compression. If neutrino cooling in 
the central region is effective enough to produce a 
strong positive temperature gradient, the helium 
core flash is ignited off center. Model mass must 
be larger than M > 1.2M Q for central helium ig- 
nition to occur (Fig. [3]). The mass Mi of the he- 
lium core when the helium core flash occurs differs 
greatly between the off-center and central ignition 
cases, as seen from the third column of Table [2] 
Central ignition occurs before neutrino cooling be- 
comes appreciable, and, hence, core masses at ig- 
nition are small: Mi < 0.4M Q . Once neutrino 
cooling becomes effective in the central region, off- 
center ignition is delayed until the core mass is 
much larger, namely, Mi > 0.5M®. 

Figure [5] shows the interior temperature as a 
function of density for models of mass 0.8 M Q 
and core mass Mi = 0.49Af Q . A constant tem- 
perature "plateau" spreads inward from the base 
of the hydrogen-burning shell. This implies that 
the thermal structure of the helium core results 
mainly from the temperature in the hydrogen 
burning shell and neutrino energy losses in cen- 
tral regions ( Fujimoto et al.l 119841 ) . Along with 
a decrease in the initial abundance of CNO ele- 
ments, the hydrogen-burning rate decreases, and, 
in order to compensate for this, the temperature in 
the hydrogen-burning shell increases, which makes 
the contribution of comprcssional heating smaller. 
Accordingly, an isothermal plateau develops in the 
region where the radiative heat transport domi- 
nates over electron condu ction for the metallic - 
ity of [Fe/H] < -5 (see iFuiimoto etail Il995|) . 
In our models, the comprcssional heating plays 
a small part to raise the maximum temperature 
in the helium zone slightly higher than the tem- 
perature in the hydrogen burning shell, as seen 
from the lowest mass model; note that for the 
massive models, the helium burning already con- 
tributes appreciably to increase the temperatures 
in the right shoulder of structure lines in this fig- 



ure. This behavior contrasts with that of model 
stars of younger population in which the compres- 
sional heating play the dominant part in determin- 
ing the maximum temperature in the helium core 
and makes it much larger than the temperature 
in the hydrogen-burning shell. Because of high 
temperature in the hydrogen burning shell, there- 
fore, Z = model stars experience central ignition 
at smaller initial masses than do model stars of 
younger populations for which the minimum ini- 
tial mass for central helium ignition is ~ 2.5M Q . 
While the density and temperature of the hydro- 
gen burning shell are in local maximum just before 
the helium ignition, hot CNO-cyclc is not still ef- 
fective. At this stage, the nuclear timescale of 13 N 
against proton capture reaction (~ 10000 sec) is 
much larger than that against /3-decay reaction 
(863 sec) and is negligible in the outcome of nei- 
ther the nucleosynthesis nor the nuclear energy 
output. 

If neutrino energy losses are sufficiently effec- 
tive before helium is ignited, the central region 
cools, and helium is ignited off-center as is the case 
in low-mass stars of younger populations. In the 
central region of the model, plasma neutrinos are 
the dominant neutrino energy-loss mechanism and 
produce a steep gradient in the rate of released en- 
ergy. On the other hand, the conduction, which is 
important in the electron-degenerate core, trans- 
port the energy towards the center where the neu- 
trino loss works. Consequently, both conductivity 
and neutrino energy-loss rates promote cooling of 
the core and, thus, delay the off-center ignition of 
helium until a larger core mass gives rise to a larger 
hydrogen-shell burning rate and a larger temper- 
ature in the hydrogen-burning shell, as seen from 
Table H 

In all of the models which experience off-center 
ignition, convection driven by helium burning ex- 
tends into the upper hydrogen-rich layers during 
the decay phase of the core helium flash, because 
of smaller ent ropy in the hyd r ogen burnin g shell, 



as shown by IFuiimoto etall (|l990l . Il995l ). The 



ingestion of hydrogen into the helium convective 
zone begins a sequence of events that leads to 
the enrichment of the surface with carbon and 



nitrogen (Hollowell et al.l 


1990; Fuiimoto et al. 


2000: Schlattl et al. 20021 


Picardi et al. 2004 


Weiss et al. 20041). Characteristics of hvdroeen 



mixing are also given in Table [2] The ingestion 
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of hydrogen occurs within a matter of days after 
the helium-burning luminosity reaches its peak. 
In the models undergoing central helium burn- 
ing, we do not find a hydrogen- mixing event; The 
core helium flash is rather weak (see in Ta- 

ble [2]), which makes it difficult for the outer edge 
of the convective region driven by heliu m burnin 
to re ach the hydrogen-containing layer (jFuiimot 
19771 ). For the models of M = 1.2M , we fol- 
low the evolution through the thermally pulsing 
AGB phase to find that the Hc-FDDM is triggered 
during the helium shell flashes. 



4. Discussion and Comparison with Other 
Works 

In this section, we compare the models in the 
literature with our models adopting similar input 
physics to see which cause the dominant effect on 
the differences in the evolution and to confirm the 
correctness of numerical computations. In addi- 
tion, we compare the models without non-resonant 
effect for a-capturc reactions to see the influence 
of uncertainty in 3a reaction rates on the stellar 
structure at Z = 0. 

4.1. Comparison with 0.8M Q Models 

A model comparable w i th ou r O.8M model is 
Model 1 of lPicardi etaD (|2004L P04), of compo- 
sition Y = 0.23 and Z = 0. In their computa- 
tions, release 4.98 of FRANEC was used and time- 
dependent convective mixing was calculated; neu- 
trino energy-loss rates by plasma-neutrino emis- 
sion were modified, with consequences being re- 
ported as minimal. Reaction rates and conductive 
opacities are, respectively, common with our mod- 
els. The neutrino energy-loss rates are common 
with our model for photo- and pair-n eutrino pro 



cesses , but they use energy-loss rate oflDicus et al 



1976 ) for bremsstrahlung and of iBeaudet et al 



19671 ) for recombination processes. For plasma 



neutrino energy losses, they adopt an energy- 
loss rate (jEsposito et al.l 120031 ) which differs only 
slightly from 196 in the temperature and density 
ranges relevant to the ignition of the helium flash. 

Although there are some differences in adopted 
input physics, the evolution of the helium core 
flash of the P04 model is similar to that of our 
08cf model. The helium core mass at the onset of 
the core helium flash is the same in both mod- 



els, namely, Mi = O.52M . When CNO-cycle 
reactions are the dominant contributors to the 
hydrogen-burning luminosity, so that the hydro- 
gen profile is very steep, the quantity Mh 6 defined 
by P04 as the mass of the helium core when the 
maximum hydrogen-burning luminosity is reached 
is nearly the same as the quantity M\ we have 
defined as the mass of the helium core when the 
core helium flash begins. The maximum helium- 
burning luminosities differ by less than a factor of 
two, being L Hc = 1-2 x 10 10 L Q in the P04 model 
and Lhc = 7.6 x 1O 9 L in our model 08cf. The 
mass at the outer edge of the convective shell at 
the onset of hydrogen mixing is the same in both 
cases, namely, O.5O6M . The values of Ai m i x (see 
Table [2]) and Xq, the mass fraction of carbon in 
the helium convective shell, are also comparable: 
Ai mix = 2.1 x 10 5 sec in the P04 model versus 
1.87 x 10 5 sec in model 08cf, and X c = 4.15 X 10~ 2 
in the P04 model versus X c = 4.28 x 10~ 2 in 
model 08cf. There is a large difference in model 
characteristics when the helium convective shell 
first appears; in the P04 model, the convective 
shell driven by helium burning appears at a mass 
shell M B cs = O.348M when L Ho = O.658L , 
while, in our model 08cf, Mbcs = O.3825Af when 
L Hc = 1-57 x 1O 2 L . The factor of 200 difference 
in the helium-burning luminosity when shell con- 
vection begins is probably simply a typographical 
error in P04, an interpretation reinforced by the 
fact that the time for the helium-burning luminos- 
ity to reach its maximum value is almost the same 
in the P04 model (723 yr) and in ours (715 yr). 
It takes more than ~ 3 x 10 4 yr for the helium- 
burning rate to increase by a factor of 200 in this 
range. 

The model of 0.8^ and Z = in FII00 is 



based on those of Fuiimoto et al. ( 19951 ) (here- 
after F95) and the comparable results are given in 
their Table 1. In F95, the values of Mf 18 *, Mbcs, 
and logig| x (L ) arc O.5116M , O.37O5M , and 
9.983, respectively. This core mass Mi coincides 
with our 08nac very closely despite the differ- 
ences in the input physics; F95 took into ac- 
co unt only the resonant 3a reaction reactions 



(jAustin. Trentelman fc Kashvl Il97lh . which is 
smaller by a factor of ~ 2.4 than the NACRE rate 
at the relevant temperature range (logT ~ 7.94). 
The smaller helium burning rate tends to delay 
the ignition of helium core flash. On the other 
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hand, the 183 formulae adopted here give larger 
conductivity than the Iben's fitting formulae used 
by F95 in the region of coulomb-liquid regime 
where the maximum temperature in the helium 
zone occurs (see fig. [T]), which works to delay the 
helium ignition due to the enhanced cooling of he- 
lium zone through the inward heat conduction in 
our model. These two effects compensate for each 
other, while the effect of larger conduction is man- 
ifest in inner ignition (or in smaller Mbcs) m the 
08nac model. In actuality, the 08cf model with the 
same input physics as the 08nac model except for 
the nuclear reaction rates results in a larger core 
mass than the 08nac model. This is because the 
cf88 3a reaction rate is smaller by 20% than the 
NACRE 3a reaction rate around the temperature 
relevant here. In any case, the dependence of Mi 
on the nuclear reaction rate is very small because 
of strong temperature dependence of 3a reaction 
rate. It i s also worth noting t hat the non-resonant 
reaction (jNomoto et al.lll985T ) has little to do with 
the ignition of helium core flash because of rather 
high temperatures in the helium zone logT > 7.9, 
although it takes a major part in the later phase 
of core hydrogen burning and the omission delays 
the depletion of hydrogen until higher tempera- 
ture is reached in the center (AT C ~ 0.04 and 
Ap c ~ 0.25). 

Since our calculation does not follow the burn- 
ing of mixed-in hydrogen, we cannot assess the 
results of FII00 and P04 with regard to the occur- 
rence of a hydrogen-burning flash in the middle 
of helium convective zone, the splitting into two 
convective shells, and the merging of the upper 
convective shell with the surface convective zone. 
This remains for further investigation and will be 
discussed in a separate paper (Suda, Fujimoto, & 
Iben, in preparation). 

4.2. Comparison with 1M Q Models 

In this section, we compare our results with 
those of FIH90, WOO, SOI, and SLL02 for the evo- 
lution to the beginning of the core helium flash 
of models of mass 1M Q and initial composition 
Z = 0. In all of the cited calculations, the core 
Hc-H flash occurs, although the size of the blue 
loop differs among the different works. 

We first compare our results with the model of 
FIH90. The distinctive feature of FIH90 model 
is the He-H shell flash during the hydrogen shell- 



burning, as stated in the introduction. We will 
show in the following subsection that the insta- 
bility of the hydrogen shell burning and the resul- 
tant shell flash arc solely attributable to the exclu- 
sion of non-resonant rate of 3a reaction reactions 
(jNomoto et al.l [h^H. Furthermore, the FIH90 



model ended in a significantly larger core mass at 
the ignition of helium core flash (Mi = O.528Af0). 
It is even larger as compared with the O.8M0 
model of FII00 with the same input physics except 
for the equation of state (EOS) , despite the gen- 
eral tendency of decrease for larger stellar masses 
as seen in Table. [U Since these two computations 
differ only in the equation of state (EOS) for a 
Cou lomb liquid and s olid among the input physics 
(see llben et al.lll992l for the adopted the EOS), 
the main reason for the larger core mass may be 
the larger radius of the helium core in the FIH90 
model (see their Table I in FIH90), resulting from 
the difference in the adopted Coulomb corrections 
in the EOS. A larger core radius implies a smaller 
gravity of the core, and hence, a smaller tempera- 
ture in the hydrogen-burning shell for a given core 
mass, to defer the ignition of helium core flash un- 
til a larger core mass is achieved. In actuality, in 
FIH90 model, a He-H shell flash is postponed un- 
til the core mass grows as large as Mi = 0.505M Q 
and ignited at a low density in the bottom of hy- 
drogen burning shell under a flat configuration 
(V » 4) but under non electron-degener ate condi- 
tions (see their Fig. 2 in Fuiimotol 1982 ). Because 
of a large core mass, the He-H shell flash has little 
effect on the thermal state of the inner core and 
FIH90 find off-center ignition of the helium core 
flash. 

Nex t, we compare with the work bv lWeiss et al. 
(|2000l . WOO) who use a code which differs 



from the one used bv fSchlattl et all (l200ll . S01) 
with regar d to the EOS (IStranierd I1988T ) . reac- 



tion rates ( Thielemann et al.lll987l . which include 
the non-resonant term in the 3a reaction rate), 
and the radiative opac i ty (old version of O PAL, 
Rogers fe Iglesiaslll992l; llglesias et al.lll992h . The 



conductive opacity and the neutrino energy-loss 
rates are not described in their paper. Since 
they do not follow the progress of the helium 
core flash, they do not find the Hc-FDDM event 



1 The computation of FIH90 was done at the University of 
Tokyo Observatory in 1984 when one of authors (I.I., jr.) 
visited it as JSPS fellow. 
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while helium is ignited off center. Their value of 
frpO = 6.31 Gyr is close to that of our model lOcf 
(6.54 Gyr); the small difference may be due to the 
use of different versions of OPAL opacities. As for 
the core Hc-H flash, their values for the location 
of the outer edge of the central convective zone 
Mecs = 0.11M© and for the helium-burning lu- 
minosity L^^ x = 2.57 x 10~ 7 Lq are similar to 
our values of M ECS = O.115M and L^ x = 

3.27 x 1O _7 L for our model lOcf. Element abun- 
dances at maximum nuclear burning luminosity 
are also comparable; WOO find Xi 2 = 6.50 x 10~ 12 , 
X u = 2.19 x 10~ 10 , and X 16 = 2.77 x 10~ 12 , and 
we find X 12 = 7.90 x 10~ 12 , X u = 1-20 x lO" 10 , 
and Xi 6 = 8.44 x 10~ 13 . These abundances 
are influenced slightly by the choice of time 
step. At the tip of the RGB, we can com- 
pare with their "canonical" model that element 
diffusion is not considered in their work; their 



tip 



O.497M and log(Lfv/L Q ) = 2.357 differ 



M 

slightly from our lOcf model: Af* ip = 0.5054M Q 
and log(L* ip /L ) = 2.413 (See Tables ED and [2]) . 

In the SOI calculations, evolution is followed 
from the main sequence to the TPAGB phase. 
The SOI input physics differs from ours with re- 
gard t o the EOS in regions of ele ctron degen- 
eracy (|Kippenhahn fc Weigertjll990h . the energy- 
l oss rates for photo, p air, and plasma neutrinos 
( Munakata et alj Il985f ) , and weak screening for 
nuclear reactions (they use the Salpeter formula). 
For the EOS in core regions, they adopt a sim- 
plifi ed equation of state for a dege nerate electron 
gas ( Kippenhahn fc Weigertl [T990h but they pro- 
vide no explicit statement as to their treatment of 
the ion gas. Their nuclear reaction rates and con- 
ductive opacities are the same as those which we 
have used. 

SOI give numerical results only for the helium 
core flash phase. Their value of log(L/L ) — 
2.314 at the start of the core helium flash is slightly 
smaller than our model lOcf value and their value 
of M-f 113 * = O.482M is also smaller than ours. 
However, they find A/g^g = O.151M , while our 
model lOcf gives Mg^g = O.339OA/ and none of 
our other models ignite a helium core flash with 
Mbcs smaller than O.3M , except in the case of 
central ignition (Table [2]). Because of the much 
larger mass between the base of the convective 
shell and the location of the base of the hydrogen- 



rich layer, the time required for the outer edge of 
the convective shell to reach hydrogen-rich mate- 
rial is much larger in the SOI model than in our 
models: Ai m j x = 10 yr for SOI versus Ai mix = 
10~ 3 - 10 -2 yr for all our models, irrespective of 
the input physics. We note that it takes more than 
1000 yr for convection generated in the center to 
reach the maximum extension in mass and Ai m j X 
for S01 falls in the middle of our two cases. 

In an effort to reproduce the S01 results, we 
constructed a l.OM model lOcf, using the fit- 
tin g formula for neutrino energy-loss rates given 
bv Munakata et all (|l985f ) which docs not include 
neutrino bremsstrahlung. The contribution of 
neutrino loss is rather small as compared with 
the gravitational energy release in the core, and 
yet, affects the internal structure of helium core 
significantly. In Fig. [5l we compare the struc- 
ture line at core mass Mi — O.49M in which 
the neutrino energy loss rate (L v = O.7OL ) is 
smaller than the rate of gravitational energy re- 
lease (L g = 2.3i ) and the helium burning rate 
is still small (Lho = O.35L ). We see that the 
maximum temperature shifts to the inner shell as 
much as A log p ~ 0.13. The reason for this inward 
shift is that, near the stellar center, the energy-loss 
rate due to neutrino bremsstrahlung is compara- 
ble to the energy-loss rate due to other neutrino 
processes; neglect of the neutrino brcmstrahlung 
contribution means that the cooling rate in cen- 
tral regions is reduced from what it would other- 
wise have been. Accordingly, the initiation of a 
helium-burning thermonuclear runaway occurs for 
a smaller core mass than would otherwise be the 
case. We find that, when Mi = O.4968M , helium 
is ignited at a mass point Mg^g = O.2933M , 
which is about 10% smaller than we find for our 
lOcf model when neutrino bremstrahlung is in- 
cluded. The 10% reduction we have found is, how- 
ever, far too small to account for the S01 result but 
we suspect that differences in neutrino energy-loss 
rates are in part responsible for the small value of 
■^/gCS found by SOI. Since radiative opacities, 
conductive opacities, and nuclear reaction rates 
are presumably also not responsible, differences in 
the EOS may be another source of the discrep- 
ancy. In dense stellar matter, Coulomb correc- 
tions reduce the pressure, leading to an increase 
in the density and to a reduction of core radius. 
The increase in the gravity entails the higher tem- 
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perature in the hydrogen burning shell and heats 
up the core. To explore this point more quanti- 
tatively, we examine conditions in the 10cf model 
for the same core mass as the SOI model ignites 
helium (Mi = 0.482M Q ). In Fig. [SJ we locate by 
filled circles the density and temperature of the 
hydrogen-burning shell and the density and tem- 
perature at the mass point M r pa 0.151M Q when 
Mi = 0.482M Q . We conclude that the SOI struc- 
ture curve must be fairly different from ours in the 
sense that their internal core might be kept much 
hotter despite the neutrino energy loss, presum- 
ably due to larger Coulomb corrections. 

Finally, we compare the 1M© model of SLL02 
with our model lOnac. SLL02 provide many evolu- 
tionary tracks of zero metallicity models covering a 
large mass range. They do not encounter the He- 
FDDM phenomenon and their evolutionary cal- 
culations for low mass stars extend to the AGB 
phase. Much of the input physics they adopt is the 
same as we have used to construct model "nac". 
The conductive opacities, radiative opacities, nu- 
clear reaction rates, and the mixing length param- 
eter (a = 1.5) are presumably the same. However, 
they use a different EOS and a dif ferent treatment 



of th e nuclear screening factor (jGraboske et al 
Il973h . They do not comment on the choice of 
neutrino energy-loss rates. 

With respect to the core He-H flash, the min- 
imum carbon abundance by mass for the appear- 
ance of convection at the center in their models is 
logXi2c — —11.5 for various model masses and 
this agrees well with our results for lower mass 
models. The main sequence lifetime, measured 
by the age at the turn-off point, is t^Q = 6.56 
Gyr in excellent agreement with our result of 
6.53 Gyr (see Table [TJ; this agreement is to be 
expected since both opacities and nuclear reac- 
tion rates are the same in both cases. For the 
same reason, the CN-cycle takes over as the main 
energy-production mechanism at essentially the 
same abundance of carbon at the center: X c = 
5.8 x 10~ 4 in the SLL02 model and 5.71 x 10" 4 
in ours. The maximum mass of the convcctivc 
core and the maximum helium-burning luminos- 
ity arc very similar in the two cases: Mg^g ~ 
0.095M Q and Lg| x ~ 10~ 7 £ Q in the SLL02 
model, compared with MgQg = 0.104Af© and 
^max ^ = 2.5i x 1CT 7 L Q in ours. At the max- 
imum helium-burning luminosity, element abun- 



dances are X 12 = 6.50xl0" 12 , X u = 2.19xl0- 10 , 
and X ie = 2.77 x 10~ 12 in the SLL02 model, com- 
pared with X 12 = 6.64 x 10" 12 , X u = 1.19 x 
10~ 10 , and X 16 = 8.97 x 10" 13 in our model. 

SLL02 do not find a shell He-H flash, consistent 
with our results. Very similar results arc obtained 
for stellar luminosity and the mass of the helium 
core at the RGB tip; they find logL^P = 2.357 
and Mj tlp = O.497M , compared with logL^P = 

2.372 and Mj lp = 0.4922M Q for our model. Core 
helium burning begins when Mg^g = 0.31Mq for 
SLL02 and M BCS = O.332OM in model lOnac 
(see Table [2]). SLL02 actually find that hydrogen 
mixes into the convective zone driven by core he- 
lium flash, but, because the hydrogen abundance 
that appears in the zone is so small (X < 10 -8 ), 
the luminosity due to hydrogen burning is rela- 
tively small (Ljj < 10 4 L©), so they neglect the 
effects. As mentioned earlier, their result may be 
an artifact occasioned by their assumption that 
convective mixing is instantaneous. 

For completeness, we discuss additional proper- 
ties of our 1M© model lOnac and, when possible, 
compare with properties of the other models. At 
the turn-off point, our model has Mi = 0.2032M Q 
and the abundance of hydrogen at the center is 
finite with the value X c = 1.57 x 10~ 2 . Af- 
ter 553 Myr of evolution beyond the turnoff 
point, a convective core is formed, driven by 
the Hc-H core flash. Still 3.06 Myr later, the 
hydrogen-burning luminosity reaches a maximum 
of ig ax = 24.6L©, with L pp = 18. 3L© and 
Lcn = 6.29L© being, respectively, the contri- 
butions of the pp-chain reactions and of the CN- 
cycle reactions. The blue loop in the H-R dia- 
gram is characterized by 3.795 < logT e g (K) < 
3.863, and 1.27 < logL s < 1.33. During the flash, 
the convective core grows to a maximum mass of 
MgQg = 0.104M©, which is between the masses 
found by WOO and SLL02. The maximum carbon 
abundance is achieved at nearly the same time 
that the convective shell achieves its maximum 
mass and has the value A™ ax = 7.146 x 10~ 12 . 
As evident from Fig. [3j the inner CN-burning 
shell passes through the site of this central con- 
vective zone when, at the center, logp c = 5.017 
and logT c = 7.801. A temperature plateau de- 
velops in the outer core, as is characteristic of 
all zero metallicity models, and as the core mass 
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continue s to increase, a temp erature inversion is 
formed ( Fuiimoto et al. 198J). Thereafter, en- 
ergy transport into the center plays a crucial role 
in determining the maximum temperature in the 
core. Our model center evolves into a region of 
strong degeneracy, reaching the maximum cen- 
tral density logp™ ax = 6.074 at logT c = 7.840. 
When the helium flash ignites off-center, the 
core mass is Mf 13 -* = O.5O28Af (see Table Q} 
and the luminosity is L max = 251. 6L Q . At 
maximum helium luminosity, the central abun- 
dances of CNO elements are X\ 2 = 1.18 x 10~ 5 , 
X u = 5.35 xlO" 10 , andXie = 9.21 xl0~ 8 . In our 
model, T c = 5.622 x 10 7 K and p c = 8.498 x 10 5 g 
cm -3 at this time. 

4.3. Influence of 3a reaction rates 

In this subsection, we discuss the relevance of 
characteristics of 3a reaction rates to the evolu- 
tion of zero-metallicity stars. Recently, the prop- 
erties of resonances other than the so-called Hoylc 
resonance in 12 C and their effects on 3a reaction 
rates have been discussed experimentally and the- 
oretically sin ce the compilation of NACRE (see 
for example. Iltoh et all 12004 iFvnbo et all l2005t 
iKurokawa fc KateTl 2005 ) . On the other hand, the 
non-resonant term of 3a reaction rate may also be 
subject to uncertainty in the cross section because 
the nuclear theory hardly determines the behav- 
ior of three body interactions (K. Kato 2007, pri- 
vate communication). Considering these improve- 
ment in the field of nuclear physics, it is impor- 
tant to examine the effects of resonant and non- 
resonant terms of 3a reactions, separately. Indeed, 
for Z = models, the non-resonant term of 3a re- 
action rate can change the evolutionary behavior 
drastically at the hydrogen shell burning phase. 

As a test for exploring the contribution of non- 
resonant term, we compute the models of 1.0 and 
1.1 by adopting the nuclear reaction rates of 



Fowler et al.l (|l975h (hereafter FCZ75) with the 



same other input physics as the models in this 
work. The main difference of FCZ75 rates from 
the NACRE rates is the consideration of non- 
resonant term in 3a reaction. Since FCZ75 have 
not yet taken into account the effect, the 3a re- 
action rate drops by far more rapidly than the 
NACRE rates below logT < 7.89. Other cross 



sectio ns may differ by within factor of 2 or 3 (jSuda 
20031 ) and do not affect the qualitative results. 



Figures [6] shows the evolutionary tracks and in- 
terior structures in the density-temperature plane 
for these models with the FCZ75 rates (refereed as 
models lOfcz and llfcz, respectively, in the follow- 
ing) and compare them with those of our models 
with the NACRE rates. The effect of different nu- 
clear reaction rates is apparent before the deple- 
tion of hydrogen in the center. Both models with 
and without the non-resonant term experience the 
core Hc-H flash, the FCZ models postpone it un- 
til higher central temperature than the NACRE 
models; in the latter models, it is ignited at the 
low temperatures where the non-resonant term 
is effective (logT < 7.89). Because of stronger 
electron degeneracy, the FCZ75 models undergo 
much stronger flashes (Lhc = 1-06 x 10 4 and 78L Q 
for models lOfcz and llfcz, respectively), as seen 
from larger loops in this figure, than our mod- 
els lOnac and llnuc (Lhc = 24 and 29L Q , re- 
spectively), which entails larger extension of flash 
convection; for models lOfcz and llfcz, the flash 
convection reaches to the shells of M conv = 0.181 
and 0.149M Q at the maximum extension, respec- 
tively, about twice as large as compared to our 
models lOnac and llnac, with the CN abundance 
of A CN = 6.0 x 10~ 10 and 2.6 x lO" 10 , respectively. 

There is also a remarkable difference in the evo- 
lutionary tracks during the hydrogen-shell burning 
between the models with and without the non- 
resonant terms. In the hydrogen burning-shell, the 
CN-burning with carbon produced by 3a reaction 
contributes considerably to the total energy gener- 
ation though the contribution to the total energy is 
still smaller than that of p-p chain reactions for the 
stars of mass 1.2 > M/Mq > 0.9, because of lower 
entropy in the hydrogen burning shell as compared 
with the stars of younger populations. In the mod- 
els without the non-resonant term, in particular, 
there form two local maxima in the energy genera- 
tion rates, a narrow one associated with CN burn- 
ing and localized at the base of small hydrogen 
abundance (X < 0.01), and a broad one associ- 
ated with pp-chain reactions in the middle of large 
hydrogen abundance (nearly half of the surface 
hydrogen abundance). When the hydrogen burn- 
ing shell passes across the inner sphere, occupied 
by the convection during the core Hc-H flash, a 
thermonuclear runaway is triggered at the base of 
hydrogen-burning shell under electron degeneracy, 
which is a similar situation as the ignition of a core 
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Hc-H flash in the center; in addition, at the base of 
hydrogen burning shell, the pressure scale height is 
much smaller than the radial distance to the cen- 
ter, i.e., V = r/\dr/dhxP\ = GM r p/rP > 1, and 
the flat c onfiguration also contribute s to the in- 
stability ( Sugimoto fc Fuiimotolll978h . It occurs 
in the shell of mass M r = 0.212M Q with Mi = 
0.340M© for mode lOfcz and M r = O.17OM with 
Mi = 0.352M o for the first shell flash of model 
llfcz. 

On the other hand, the models with the non- 
resonant term stably burn hydrogen, although 
electrons are degenerate ~ 8) at the base of the 
hydrogen burning shell. This different behavior 
stems mainly from the difference in the tempera- 
ture dependence of the resonant and non-resonant 
terms, rather than from the difference in the burn- 
ing rate itself, as can be seen from the analy- 
sis of thin shell burning bv lSchwarzschild fc Harml 
(|1965l . see also Fujimoto 1982). At temperatures 
of logT < 7.9, the 3a reaction rates by FCZ75 
and NACRE give a large difference in their tem- 
perature dependences when compared at the same 
burning rates. For less massive FCZ75 models of 
M < 0.9M Q , the entropy is smaller and hydrogen 
is burnt before carbon production becomes appre- 
ciable. For FCZ75 models of mass M > 1.2M Q , 
higher temperature as well as weaker electron de- 
generacy tend to stabilize hydrogen shell burning; 
in addition, as the contribution of CN burning 
with products of 3a reaction reaction as cataly- 
sis overweighs the p-p chain reactions, the pres- 
sure scalehcight grows comparable to the radial 
distance to the center, whicl also stabilize the 
shell burning by making the heat capacit y negative 



with hydrostatic readjustment (e.g., see iFuiimoto 
19821) . It is noted that the temperature depen- 
dence of resonant reaction rate decreases with in- 
crease in the temperature, which also stabilizes 
the hydrogen shell-burning in combination with 
the reduction in electron degeneracy when the core 
grows more massive than Mi > O.4M , as is the 
case for lOfcz model. 

During the shell flash, the CN-cycle burning 
shell expands, and the hydrogen exhausted core 
also expands due to the reduction in the weight 
of overlying layers; the central temperature and 
density decrease almost adiabatically (Label "A" 
in Fig. 0. During the decay phase of the first 
shell flash, after the flash-driven convcctivc zone 



disappears, the core is heated by the flow of en- 
ergy from the burning shell. For model llfcz, the 
second shell flash is ignited at the shell of mass 
M r = 0.287M o with M x = O.374M , when the 
hydrogen shell-burning passes across the shells, in- 
corporated into the convection and enriched in CN 
elements during the first shell flash. This flash 
grows so strong as to drive the convection deep 
into the hydrogen-rich envelope up to the shell 
of Mi = 0.404M Q (Label "B"), which expedites 
the growth of helium core to cause the increase 
in central density. Since the matter in the con- 
vective zone driven by the shell flash is enriched 
in CNO elements, the CN-cycle reactions domi- 
nate the energy generation rate while the burning 
shell traverse in the site of convective zone. Dur- 
ing this stable burning phase, the growth rate of 
core is large because of small hydrogen abundance 
therein. Accordingly, concomitant rapid compres- 
sion of helium core increases the central tempera- 
ture. This enhanced growth of helium core after 
the shell flashes leads directly to the ignition at 
the center of a helium core flash (Label "C" ) with 
a small core masses. This is the case if the con- 
vective zone is sufficiently large, as in the case for 
the second shell flashes of llfcz model, while it is 
barely missed in lOfcz model. 

The most important point by this test is that 
the contribution of the non-resonant term is dis- 
cernible in the circumstance of Z — 0, although 
it may be difficult to detect differences by the 
observations. There is a difference in the mass 
range that the off-center helium core flashes oc- 
cur. If the non-resonant term is included in the 3a 
reaction rate, the central helium burning occurs 
at > 1.2M . Otherwise, it occurs at > I.OMq 
if we use the smaller conductivity by Ib en's ap- 



ty bylbe 

proxima t es to iHubbard fc Lampel 1 1969^ and by 

i. ( 



Canutol (I1970f ) and the neutrino loss rates by 



Beaudet et al.l (|1967n as in FII00, both slightly 
smaller than ours. Since the lifetime of 1.0M Q 
model is ~ 7 Gyr (Table [l}, such stars cannot be 
seen in the present halo if they were born. Only 
clue to those objects will be binary mass trans- 
fer between giants of mass > I.OMq and dwarfs 
of mass < O.8M . Since FII00 predict the abun- 
dance ratio of C/N - 1 for He-FDDM at RGB, 
while C/N > 5 for He-FDDM at AGB, it will be 
crucial for the constraint on the estimate of nu- 
clear reaction rate to determine the abundance ra- 
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tio of C/N and the mass of the primary for Z = 0. 
5. Conclusions 

We have explored the evolution of low-mass, 
zero-metallicity stars with the most recent input 
physics. The mass range is 0.8M Q < M < 1.2M 
in step of 0.1 M Q and the initial composition is 
X = 0.767, Y = 0.233, and Z = 0. Calculations 
extend from the zero-age main sequence to the be- 
ginning of hydrogen mixing into the helium con- 
vective region on the RGB or at the start of the 
TPAGB phase, depending on the characteristics 
of the helium core flash. 

1. The emergence of CN-cycle reactions as im- 
portant contributors to nuclear energy pro- 
duction occurs during the core hydrogen- 
burning phase in models of mass M > 
0.9M Q in consequence of the formation 
of carbon by the 3a reaction. This phe- 
nomenon is independent of the adopted 
physics and its importance is a function only 
of the initial mass and metallicity. 

2. The models of M < I.IMq undergo an off- 
center helium flash and hydrogen-mixing 
into helium flash-convection, leading to 
helium-flash driven deep mixing at the tip 
of red giant branch. On the other hand, 
the models of mass M > 1.2M Q ignite 
the core helium flash at the center and 
postpone the He-FDDM until the helium 
shell flashes occur during the early phase of 
thermal pulsation at the asymptotic giant 
branch. For models of mass 0.8Mq and 
0.9M Q , our results coinci de qualitatively 



with those first found by IFuiimoto et al 



( 2000h . 



For 



(|l99dh and IFuiimoto et al 
I.OMq < M < I.IMq, whether or not a he- 
lium core flash is ignited off the center and 
hydrogen is mixed inward into the convec- 
tive zone driven by it depend on the adopted 
nuclear reaction rates. 



We have also compared our results with those 
of other investigations. Our models made with the 
most up-t o-date input physics agree well w i th the 



mo dels oflPicardi et alJlEool IWeiss et all (|2000h 



and Siess. Livio fc Lattanziol ( 2002h during evolu 



Siess. Livio fc Lattanzio (|2002l) . although we do 
not follow evolution after the mixing of hydro- 
gen into the helium flash driven convective zone 
at the beginning of the helium core flash. To 
check the behavior of the He-FDDM event, we 
need to treat mixing with a time-dependent al- 
gorithm. On the other hand, even after adopt- 
ing the same radiative opacities, conductive opac- 
ities, nuclea r reaction rate s , and neutrino energy- 
loss rates as ISchlattl et al.l (|200l[ ) , we are not able 
to obtain the inner ignition at the onset of core 
helium burning which they find; we suspect that 
the discrepancy may be due to differences in the 
EOS for the Coulomb corrections in the liquid and 
solid states and due to the neglect in their work 
of neutrino energy losses associated with neutrino 
bremsstrahlung. 

By treating the resonant and non-resonant 
rates of 3a reactions separately, we demonstrate 
that the non-resonant term plays a critical role 
in the low-mass, zero-metal stars. The neglect of 
non-resonant term causes the lower border in mass 
of central helium burning, i.e., > I.OMq. This ex- 
plai ns the discrepan c y of th e results between ours 



and IFuiimoto et al 



(|2000h . which stems mainly 
from the difference in the temperature depen- 
dence rather than in the energy generation rates 
themselves. We first point out the possibility of 
discerning the effect of non-resonant term of 3a 
reaction from the evolution of stars other than at 
low temp erature regime in th e accreting degener- 
ate stars (jNomoto et al.lll985l ). It is important to 
precisely determine the abundances and the prop- 
erties of extremely metal-poor stars to constrain 
the nuclear reaction rates. 

We are grateful to I. Iben Jr. for improving 
and revising our manuscript. We wish to thank A. 
Ohnishi and K. Kato for valuable comments on un- 
certainties on nuclear reaction rates. This work is 
part of a PhD. thesis constructed at Hokkaido Uni- 
versity and is in part supported by a Grant-in- Aid 
for Science Research from the Japanese Society for 
the Promotion of Science (15204010, 18104003). 



tion on the main sequence and the RGB. In par- 
ticular, we obtain nearly the same results as do 
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A. Conductive Opacities 

The most up-to-date conductive opacities are those of Itoh et al. ( 19831 hereinafter 183). We consider an 
ion mixture consisting of n species of nuclei. The conductive opacity n c for temperature T§ (units of 10 8 K) 
and density p 6 (units of 10 6 g cm -3 ) is taken from cq. (7) in 183, 



1.280 x 1(T 3 I ^XiAi (S) i 



1 + 1.018 



2/3 



Ai 



2/3 

Pe 



cm g | 

Pe/ 



(Al) 



where () denotes the average over the nuclear species and Zi, Xi, and Ai are the atomic number, the mass 
fraction, and the atomic mass number of i th nucleus. The average value is taken as 



{S ~ l}i i + i,oi8(EW) 2/3 pf ( Yr 



(A2) 



The quantities (S-i) i and (S+-\ ) - are c alculated with eqs. (8) and (9) in 183 and with the parameters 
according to eq.(19) of lltoh et all (|2004h : 



a e ksT 



0.2275-' 



Z 



5/3 



1/3 



PeY, 



3 = 1 



X 3 Z 3 
A, 



0.45641 lnTi - 1.31636 



-1/3 



1.388 x 10" 2 ( -^ X 3Pe 



(A3) 
(A4) 

(A5) 



J =1 



where Ti is the Coulomb coupling constant for the i th nucleus, a e the electron-sphere radius defined as 
a e = (3/47m e ) 1 / 3 with the electron number density n e , and r s the electron density parameter. 

The expressions for various mixtures are to be be considered as first approximations, compared with the 
formulae in 183 which are for pure compositions and are accurate solutions. In particular, the approximation 
for a mixture of elements with very different Z's (e.g., 1 H and 56 Fe) is not very accurate; however, the exact 
solution for such mixtures is not presently attainable. For mixtures of elements of comparable Z's (e.g., 12 C 
and 16 0), the approximation is fairly accurate. 

The 183 results are strictly applicable only for T < Tp, where Tp is the Ferm i temperature defined by 
eq. (1) in 183. In this work, defined by 183 and defined bv llbenl (Il975h are interpolated with one 

other over the range 0.5 < T/Tp < 2.0 using a sin squared algorithm. 

Another algorithm is used to extrapolate 183 results for values of T outside the region defined by 2 < T < 
160. Both (S-i) and (£+i) are constrained by (5_i) < (<5-i) lim and (5+i) < (5 + i) lim , where the upper 
limits arc calculated by demanding that, in eq. (6) in 183, 



k 
2kp~ 



1, 



k 
2k~ P 



,0 



k 2 

""TP 



k 2 



(A6) 
(A7) 



where /ctf is the Thomas-Fermi wavenumber. The first condition assumes that interactions between ions can 
be neglected and the second condition expresses the Thomas-Fermi approximation for electron screening. 
These conditions give for the upper limits: 



(S-i) u 



In 1 



4^ 



- 1 



k 2 

44 



(A8) 
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(S + i) l 



k 2 



In 1 



The factor kj, F /{4k F ) is given by 



k 2 



44 



1 



2ll 

1 

137' 
mc 



k 2 



1 

137 



1/2 



1/3 



A- 2 



44 



(A9) 

(A10) 
(All) 
(A12) 



where a is the fine structure constant and 6 is taken from 183. This holds both for relativistic degeneracy 
(P6 > 1) and for non- relativistic degeneracy (pe < 1). 

The interpolation between 183 and Hubbard fc Lampel ( 19691 ) as approximated by llben (1975) is shown 
in Fig.[TJ In the actual computations, the upper limits lA"8l and lAUl are infrequently invoked. 



B. Neutrino Losses By Bremsstrahlung 



We describe here the est imate of the neut rino energy-loss rate due to bremsstrahlung for the case of an 
ion mixture as discussed bv lltoh et al. ( 19961 . hereinafter 196). For any nuclear species, gas, liquid, and solid 
states for ions are taken into consideration. The total energy-loss rate is given by 



<?br 



£^(Qgas + Qii q + Q 



'sol 



(Bl) 



where the Xi are nuclear abundances and the summation is taken over the species 4 He, 12 C, and 16 0. The 
gas state corresponds to the scheme for weakly degenerate electrons and Qg as is based on eq. (5.1) in 196: 



Z 



Qgas = 0.5738 -f-2* p (C + F gas - C_G gas ) , 



(B2) 



where 



C- = \{((%-CP A )+n(C9-C%)} 



(B3) 
(B4) 



All values of the right hand side of these equations are defined in 196. Fgas an d Ggas are independent on 
the nuclear species other than the number of electrons. -Fgas is given by eq. (5.2) in 196, but the parameter 
77 is given by 



^ = £^ p(7 . 05xl0 6 T 1.5. 



5.12 x 10 4 T S 3 1 



(B5) 



Ggas is given as follows according to eq. (5.9) in 196 
G 



'gas 



1 + 10 -9£^ p ) (a3 + a4Ts 



-2 
8 



as T 8 



EX, Zi 
— —p 



A, 



0.656 



(B6) 
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According to 196, the applicable range of Qg as is determined by T > 0.01 Tp for 4 He and T > 0.3 Tp for 
12 C and 16 O, where Tp is the Fermi temperature: 



T F = 5.9302 x 10 9 { \jl + 1.018 (jl^^-Pej -l\ [K] . (B7) 

Qj iq is calculated according to §5.2 in 196. For T < 0.3 Tp (or T < 0.01 Tp for 4 Hc), the formula for the 
bremsstrahlung pr ocess in th e liqui d state is adopted if the parameter, which considers ion-ion correlation 



effect according to lltoh et al.l (|2004l ) , 

Z? /3 / A XiZ. 



^ = 0.2275 ^-^^^1 (B8) 



satisfies F; < 180. Here, is given in cgs units as 



2liq = °' 5738 P ( C + F hq - C - G hq) ' (B9) 

where and Gjj are provided by cqs. (5.19) and (5.20), respectively, in 196. These quantities are 
functions of u and Ti via 

3 

vi = 5Z« m r^ f , (bio) 

m=0 
3 

Wi = ^/3 m Fr f - (Bll) 

m— 

For Ti > 180, we assume for the bremsstrahlung process that 

Q sol = ^lattice + ^phonon ( B12 ) 

where, according to eq. (5.29) in 196, 

Z 2 

^lattice = °- 5738 ^" T 8 ^ (C+^lattice ~ C - G lattice) ' ( B13 ) 

Qphonon = °- 5738 f T 8V(c + F phonon -C_G phonon ). (B14) 



Note that we adopt /b anc [ = 1 in QJ a ^j ce - Each factor such as F and G in the right-hand side of IB13I and 
IB 141 depends on the nuclei via the parameters: 

3 3 

«< = $Z « m r i " f Wi = (B15) 

m— m— 

3 3 



m— m— 
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Table 1: Evolutionary Characteristics of Low- 
Mass, Population III Models with Up-to-Date In- 
put Physics 
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7.51 x 10 


-3 


lOnac 


0.5028 


0.3320 


9.834 


237.1 


9.038 


846 


9.81 x 10 


-3 


llnac 


0.4929 


0.2986 


9.722 


226.9 


8.744 


953 


1.82 x 10 


-2 


12nac 


0.3840 


0.0 


3.829 


126.8 










08cf 


0.5151 


0.3717 


9.882 


248.6 


9.264 


715 


5.91 x 10 


-y 


09cf 


0.5121 


0.3629 


9.858 


245.3 


9.231 


749 


6.49 x 10 


-3 


lOcf 


0.5054 


0.3390 


9.831 


239.4 


9.026 


858 


9.99 x 10 


-3 


llcf 


0.4966 


0.3115 


9.728 


231.1 


8.776 


941 


1.69 x 10 


-2 


12cf 


0.3829 


0.0 


3.393 


124.2 











Table 2: Characteristics of Helium Core Flash of 
Model Stars with Up-to-Date Input Physics 
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log T (K) 



-0.5 



Fig. 1 . — Inter polation of conduct i ve op acities 
between 183 and Hubbard fe Lampel ( 19691 ) as ap- 



proximated by llbenl (|l975h ~ for pure He com- 
position. Conductive opacities by 183 (dotted 
lines) and by Iben's fitting formulae (dashed lines) 
are interpolated (solid lines) in the range 0.5 < 
T/Tp < 2.0, where Tp is the Fermi temperature. 
Each Line gives the conductive opacity for a con- 
stant density overthe range 10 3 - 10 7 g cm -3 . The 
shaded area surrounded by dotted line represents 
where the conductivity by 183 becomes effective 
for density greater than given by Iben's fitting 
formulae. Thin solid line shows the distribution 
of opacities in the lOnac model during the helium 
core flash function of temperature. 
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Fig. 2. — H-R diagram for models of masses in the 
range from 0.8 to 1.2M Q in steps of 0.1M Q with 
the nuclear reaction rates g iven by NACRE Of; 
thick lines] aug99 and by ICaughlan fe Fowler 
( 19881 . ; thin lines). Evolutionary tracks are 



plotted from the zero-age main-sequence to the 
hydrogen-mixing event at the tip of RGB (for 
O.8M > M > 1.1M©) and to the early thermal- 
pulsating AGB (for M = 1.2M Q ), respectively. 
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Fig. 3. — Evolutionary tracks of model centers 
in the temperature-density plane. The models are 
the same as in FigJ^l 



Fig. 4. — Variations in several physical quantities 
as a function of the hydrogen abundance at the 
center. In order from top to bottom, panels show 
the central temperature, mass fraction of the sum 
of CNO elements at the center, and the mass of the 
convective region. Each line represents the time 
variation of 0.8 (solid lines). 0.9 (dashed lines), 
1.0 (dot-dashed lines), 1.1 (dotted lines), and 1.2 
Mq (short-dashed lines) model, respectively. 
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Fig. 5. — Stellar structure in the temperature- 
density plane for various models at an evolution- 
ary stage corresponding to a helium core mass 
of Mi — 0.49Af Q . Model names are described 
in Table [1] The two filled circles show p and T 
at the hydrogen-burning shell and at a shell of 
mass M r = 0.151M Q in the 10cf model for which 
Mi = O.482A/0, the same core mass as in SOI at 
the onset of helium ignition. For comparison, the 
thin solid line describes the lOcf model a few Myr 
before helium ignition. Triangles on the line lo- 
cate the ignition point (right hand side) and the 
hydrogen-burning shell (left hand side). 



Fig. 6. — Evolutionary tracks (thick lines) 
and structure at some stages (thin lines) in the 
temperature-density plane for selected models. 
The labels on top left corner denote the mass 
and the adopted nuclear reaction rates for each 
model. Circles and triangles denote the center of 
each model and the hydrogen-burning shell where 
the hydrogen-burning rate is maximum. The solid 
straight line and the dashed straight line denote 
constant T = 0.5 and T = 0.3Tp, respectively 
where T is the Coulomb coupling constant and Tf 
is the Fermi temperature. See text for detail on 
labels "A" to "C". 
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